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Abstract 

Small  particles  combusting  under  free  molecu¬ 
lar  conditions  are  shown  to  exhibit  critical  behavior 
under  conditions  of  high  collision  efficiency.  At  a 
particular  collision  efficiency  the  particle  behavior 
transitions  from  non-critical  to  critical  (or  runaway) 
behavior.  Non-critical  behavior  is  defined  as  a  finite 
particle  temperature  at  burnout,  while  critical 
behavior  implies  an  infinitely  increasing  particle 
temperature  at  burnout.  The  critical  collision  effi¬ 
ciency  is  derived  from  the  governing  equations  for 
free  molecular  combustion  of  small  particles.  The 
critical  collision  efficiency  is  shown  to  be  a  function 
of  the  specific  heats  of  the  gas  and  the  particle, 
and  is  independent  of  the  particle  size.  Analytic 
and  numerical  results  are  compared.  Particle  radi¬ 
ation  is  included  in  the  numerical  results.  The  inclu¬ 
sion  of  radiation  does  not  effect  the  critical  collision 
efficiency,  although  radiation  does  decrease  the 
rate  of  particle  temperature  increase. 

Nomenclature 

A  Particle  surface  area 

C  Particle  specific  heat 

Cp  Gas-specific  heat  at  constant  pressure 

Gas-specific  heat  at  constant  volume 

Cy  Gas-specific  heat  at  constant  volume 

under  effusive  conditions 

D  Particle  diameter 

E  Particle  internal  energy 

e  Particle  specific  internal  energy 

f  Collision  efficiency 


h°{0)  Enthalpy  of  formation  at  0  K 

k  Boltzmann  constant 

Spectrally  averaged  emission  coefficient 

m  Gas  molecular  mass 

m  Mass  of  carbon  atoms  removed  from  parti¬ 
cle  per  reactive  collision 

M\^  Gas  molecular  weight 

n  Gas  number  density 

Na  Avogadro’s  number 

PU  Ambient  pressure 

q„r  Energy  transferred  to  the  particle  by  a  non¬ 
reactive  collision 

Energy  transferred  to  the  particle  by  a 
reactive  collision 

9t  Specific  gas  constant 

r  Particle  radius 

/q  Initial  particle  radius 

t  Time 

t*  Particle  burnout  time 

T  Particle  temperature 

Tq  Reference  temperature  or  initial  particle 
temperature 

7^  Ambient  temperature 

"f  Particle  burnout  temperature 

u  Gas  thermodynamic  internal  energy 

V  Particle  vblume 

v  Average  molecular  velocity 

X  Stoichiometric  coefficients 


*  The  research  reported  herein  was  performed  by  the  Arnold  Engineering  Development  Center  (AEDC),  Air  Force  Materiel 
Command.  Work  and  analysis  for  this  research  were  performed  by  personnel  of  Sverdrup  Technology,  Inc.,  AEDC  Group,  technical 
services  contractor  for  AEDC.  Further  reproduction  is  authorized  to  satisfy  needs  of  the  U.  S.  Government. 
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Z  Collision  frequency 


E=pV{t)eiir)  (1) 


a  Particle  surface  thermal  accommodation 
coefficient 

Y  Gas  ratio  of  specific  heats 

e  Particle  emissivity 

X  Gas  mean  free  path 

p  Gas  kinematic  viscosity 

p  Particle  material  density 

Ambient  gas  density 
o  Stefan-Boltzmann  constant 

(0  Specific  oxidation  rate 

Introduction 

The  term  “critical  behavior”  is  used  to  describe 
many  different  aspects  of  combustion,  e.g.,  ignition 
vs.  non-ignition,  explosion  vs.  non-explosion,  vigor¬ 
ous  vs.  weak  combustion.  In  general,  criticality 
indicates  an  abrupt  transition  from  a  stable  or  well- 
behaved  configuration  to  a  completely  different, 
potentially  unstable  configuration.  The  critical  com¬ 
bustion  behavior  of  spherical  carbon  particles  has 
been  examined  previously  by  Trevino,  Higuera, 
and  Linan.''  Additionally,  Zhang^  analyzed  bifurca¬ 
tion  behavior  in  a  homogeneous-heterogeneous 
system.  In  each  of  these  instances,  the  particle 
combustion  was  assumed  to  occur  under  contin¬ 
uum  conditions.  Hiers®’'^  has  shown  that  the  com¬ 
bustion  of  soot  particles  occurs  under  substantially 
rarefied  or  free  molecular  conditions  for  most  prac¬ 
tical  applications.  The  objective  of  this  paper  is  to 
identify  a  regime  of  critical  behavior  for  particles 
combusting  under  free  molecular  conditions  -  indi¬ 
cated  by  an  infinitely  increasing  particle  tempera¬ 
ture  at  burnout. 

Free  Molecular  Combustion  Model 


where  p  is  the  material  density  (assumed  con¬ 
stant),  V  is  the  particle  volume,  and  e  is  the 
particle-specific  internal  energy  (energy  per  unit 
mass).  Since  particle  volume  is  simply  4/3k.i^,  then 
Zis  a  function  of  both  particle  radius,  r,  and  particle 
temperature,  T.  Since  both  V  and  T  (and  thus  e) 
are  functions  of  time,  then  the  rate  of  change  of  Eis 

dt  ~  dedt  dVdt  ^  ’ 


The  partial  derivatives  are  obviously 


dE 

de 


p/ 


and 


dE 

dV 


=  pe 


(3) 

(4) 


since  we  are  treating  the  particle  material  density 
as  a  constant.  Additionally,  we  can  apply  the  chain 
rule  to  yield 


where 


W 

dr 


dt  dr  dt 

(5) 

5^1 

II 

II 

(6) 

where  A  is  the  particle  surface  area.  Making  these 
substitutions  into  Eq.  (2)  yields 


dt 


.,de  ^.dr 


(7) 


We  can  write  the  internal  energy  as 
T 

e  =  J  Cdf  +  e{  (8) 


Complete  details  of  the  free  molecular  combus¬ 
tion  model  are  given  in  Refs.  3  and  4.  A  brief  deri¬ 
vation  is  given  here.  Each  term  is  defined  in  the 
Nomenclature  above.  All  thermal  properties  of  the 
particle  and  the  surrounding  gas  are  assumed  to 
be  temperature-dependent  in  this  derivation.  Addi¬ 
tionally,  the  particle  is  treated  as  a  lumped,  one¬ 
dimensional  spherical  system.  If  we  denote  the 
particle  thermal  energy  as  E,  then  Eis  given  by 


where  C  is  the  particle  specific  heat  (a  function  of 
temperature)  and  Tq  is  some  reference  tempera¬ 
ture.  If  we  use  the  thermodynamic  properties  of 
graphite  (a  reference  element)  for  soot,  then  e(To) 
is  defined  to  be  zero.  Throughout  this  paper,  the 
reference  temperature  (Tq)  is  assumed  to  be  0  K. 
Applying  Leibnitz’s  Rule  yields 
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7(0 

=  cinoi|r(o-cr4r,+  J  ^ 

(9) 

{CCT)]df 

but 

II 

(10) 

and 

IfiCCm  =  0 

(11) 

since  C  is  not  an  explicit  function  of  time.  There¬ 
fore,  Eq.  (9)  simplifies  to 

II 

(12) 

Substituting  Eq.  (12)  into  Eq.  (7)  yields 

dE  i/^dT  ^.dr 

Tt  = 

(13) 

Furthermore,  since 

(14) 

then  Eq.  (13)  will  simplify  to 

dt  ~  dt  ^  dt\ 

(15) 

If  we  assume  we  are  given  the  specific  oxida¬ 
tion  rate,  ©(g/cm^/sec),  as  a  function  of  tempera¬ 
ture  and  pressure,  we  can  find  the  rate  of  change 
of  the  particie  radius  since 

Q- 

II 

Q. 

II 

3 

1 

(16) 

SI  CL 

1 

II 

o 

(17) 

Substituting  Eq.  (17)  into  Eq.  (15)  yields 


This  is  the  rate  of  change  of  the  energy  content  of 
the  particle. 


The  free  molecular  model  is  one  of  collision¬ 
counting.  A  collision  which  results  in  an  exothermic 
reaction  adds  energy  to  the  particle.  An  endother¬ 
mic  reaction  removes  energy.  Non-reactive  colli¬ 
sions  will  either  add  or  subtract  energy  from  the 
particle:  adding  if  the  gas  is  hotter  than  the  particle, 
subtracting  if  the  particle  is  hotter  than  the  gas. 
Also,  if  a  reaction  occurs,  mass  must  be  removed 
from  the  particle. 

If  the  characteristic  body  dimension  is  small 
compared  to  the  gas  mean  free  path,  the  body  is  in 
free  molecular  flow,  and  heat  transfer  will  occur 
only  by  direct  collisions  with  the  ambient  gas  mole¬ 
cules  (ignoring  radiation  for  the  moment).  Upon 
collision  with  the  particle  surface,  the  gas  mole¬ 
cules  will  exchange  energy  with  the  surface,  either 
by  transferring  thermodynamic  internal  energy  and 
translational  energy,  or  by  chemically  reacting  with 
the  particle  surface.  The  term  "thermal  accommo¬ 
dation"  refers  to  the  effectiveness  of  the  energy 
exchange  between  the  gas  and  the  particle  during 
a  collision.  Perfect  thermal  accommodation  is 
assumed  throughout  this  work,  which  means  that 
the  energy  modes  of  a  gas  molecule  come  to  ther¬ 
mal  equilibrium  at  the  particle  temperature.  This  is 
reasonable  for  low  thermal  energy  collisions  with 
rough  surfaces  such  as  those  considered  here.® 
We  can  express  the  oxidation  rate  (co)  as  a  fraction 
of  the  total  number  of  collisions;  i.e., 

to  =  Zfm  (19) 

where  Z  is  the  total  collision  frequency,  Zf  is  the 
frequency  of  reactive  collisions,  f  is  termed  the  col¬ 
lision  efficiency,  and  fn  Is  the  mass  removed  per 
reactive  collision.  Assuming  no  relative  motion 
between  the  particle  and  the  gas  (dynamic  equilib¬ 
rium),  the  collision  rate  Zis  given  by  Ref.  6 


where  ^is  Boltzmann’s  constant,  v  is  the  average 
molecular  velocity,  n  is  the  number  density,  m  is 
molecular  mass,  and  and  are  the  pressure 
and  temperature  of  the  ambient  gas.  The  energy 
flux  at  the  particle  surface  due  to  collisions  with  gas 
molecules  can  therefore  be  given  as 


0^ 

dt  collisions 


[Zfq,  +  Z{1-f)q„,]A  (21) 
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where  q/\s  the  energy  transferred  to  the  particle  by 
a  reactive  collision  and  q„r  is  the  energy  trans¬ 
ferred  by  a  non-reactive  collision.  The  internal 
energy  {0)  of  a  gas  molecule  has  contributions 
from  translational  energy,  internal  structure  (e.g., 
vibrational  and  rotational  energy)  and  chemical 
potential  energy  (bond  energy).  We  can  write  the 
internal  energy  of  the  gas  molecules  as 

^  ~  ^translational'^  ^ Internal"^  ^chemical 


Since  the  particle  is  smaller  than  the  mean  free 
path,  the  molecular  flux  to  the  particle  surface  is 
effusive.  Therefore,  the  average  translational 
energy  of  the  gas  molecules  incident  on  the  parti¬ 
cle  surface  is  2^r/Also,  since  the  reference  tem¬ 
perature  is  chosen  as  0  K,  the  chemical  potential 
energy  is  simply  the  enthalpy  of  formation  at  0  K, 
or  h°f{  0) .  Therefore,  the  internal  energy  of  a  gas 
molecular  at  any  temperature  7"is  shown  in  Ref.  8 

T 

u{T)  =  2kT+\  interna!^  ^  (23) 

0 


The  change  in  the  internal  energy  during  a  sin- 
gie  gas-surface  coiiision  is  given  by 

^ products~  ^ reactants  (2^) 

or,  using  Eq.  (25), 


Ai/=  X  Xi 

f  products 


Ti 

J"  ^v,  Internal^ r h ^{O) 


(25) 


- 

}  reactants 


2k Tj  +  j  interna!^ h ^{O) 


q^  =  -Ai/^  (27) 

The  JANAF  Thermochemical  Tables®  and 
associated  curve  fits^®  are  used  for  all  specific 
heats  and  heats  of  formation  appearing  in  Eq.  (25). 


Equation  (21)  represents  the  flux  of  energy  at 
the  particle  surface  due  to  collisions  with  ambient 
gas  molecules.  The  particle  will  also  exchange 
energy  with  the  environment  by  radiation.  The  rate 
of  radiant  energy  transfer  is  given  by 


dt  radiation 


(28) 


where  a  is  the  Stefan-Boltzmann  constant.  The 
treatment  of  small  particle  radiation  is  explained  in 
more  detail  in  Appendix  D  of  Ref.  3.  Briefly,  small 
particles  are  mass  (or  volume)  radiators. The 
emissivity  (e)  is  thus  directly  proportional  to  the 
particle  radius.  Therefore, 

^  =  rkJT'^-T^']A  (29) 

at  radiation  ^  oo  j 

where  (the  spectrally  averaged  emission  coeffi¬ 
cient)  is  a  function  of  particle  temperature,  and  is 
derived  from  the  Lee  and  Tien^^  soot  optical  prop¬ 
erties.  The  total  energy  flux  at  the  particle  surface 
is  therefore  the  sum  of  Eqs.  (21)  and  (28),  or 


Substituting  Eq.  (19)  into  Eq.  (18)  yields 


dt  fluxes 


(31) 


where  X  denotes  the  stoichiometric  coefficients  of 
the  reaction  under  consideration.  Equation  (25) 
can  be  used  to  find  the  change  in  internal  energy 
for  either  reactive  or  non-reactive  collisions.  By  the 
First  Law  of  Thermodynamics,  the  energy  trans¬ 
ferred  to  the  particle  during  the  collision  is  simply 
the  negative  of  the  change  in  internal  energy  of  the 
collision.  Therefore,  for  a  non-reactive  collision 

dnr  ~  ~^^nr  (2®) 

and  for  a  reactive  collision 


Equating  Eqs.  (30)  and  (31)  and  solving  for  the 
rate  of  change  of  temperature  yields 


dT 

dt 


3Z 

prC 


[f{me+  q;)  +  (1-f)q„;\ 


(32) 


Note  that  the  particle  radius,  r,  appears  in  Eq.  (32) 
as  an  unknown  function  of  time.  However,  substi¬ 
tuting  Eq.  (19)  into  Eq.  (17)  yields  an  ordinary  dif¬ 
ferential  equation  for  the  particle  radius 
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%  =  ^ 
at  p 

Since  f  is  also  a  function  of  particle  tempera¬ 
ture,  Eqs.  (32)  and  (33)  form  a  coupled  set  of  non¬ 
linear  ordinary  differential  equations  for  the  particie 
temperature  and  radius  as  functions  of  time.  These 
equations  are  statements  of  energy  conservation 
and  mass  conservation,  respectively.  In  later  sec¬ 
tions  these  equations  will  be  integrated  in  time 
using  a  modified  implicit  Euler  method  described  in 
Appendix  E  of  Ref.  3. 


where 


and 


nt)  =  -p  +  (ro+ 


[fqr+{1-f)CvTJ 

[fmC-i1-f)Cv] 


5  = 


i1-f)Cv 
fmC  . 


and  7"^7is  the  initial  particle  temperature. 


(38) 


(39) 


(40) 


Critical  Combustion  Behavior  Analysis 

Equations  (32)  and  (33)  comprise  a  first-order, 
nonlinear  set  of  ordinary  differential  equations  for 
the  particle  radius  and  temperature  as  functions  of 
time.  The  nonlinearity  is  due  both  to  the  radiation 
term  ( 7^  and  the  noniinear  dependence  of  the  var¬ 
ious  coefficients  (collision  efficiency,  specific  heat, 
etc.)  on  the  particle  temperature.  The  equations 
can  be  linearized  by  neglecting  radiation  and 
assuming  that  all  thermal  properties  are  constant.  If 
the  thermal  properties  are  assumed  constant,  then 

e=  CT  (34) 


This  analytic  solution  was  used  to  validate  the 
numerical  method.  For  these  calculations,  the  con¬ 
stant  coefficients  are  evaluated  at  the  following 
conditions:  pure  moiecular  oxygen  at  a  pressure  of 
0.10  atm,  temperature  of  2500  K,  using  oxidation 
rates  from  Roth,  et  al.^®  Figure  1  presents  a  com¬ 
parison  of  the  analytical  solution  and  the  numerical 
solution  [solving  Eqs.  (32)  and  (33)]  using  the  mod¬ 
ified  implicit  Euler  scheme  described  in  Appendix  E 
of  Ref.  3  for  the  temperature  rise  of  particles  of  var¬ 
ious  size  oxidizing  under  the  conditions  given 
above.  The  numerical  solution  has  been  altered  to 
neglect  radiation  and  to  use  the  same  constant 
properties  as  the  analytic  solution. 


and 

Qnr  ~ 

=  Cy{T„-T) 

from  which  the  definition  (the  gas-specific 
heat  at  constant  volume  under  effusive  conditions) 
is  obvious.  Making  these  substitutions  in  Eq.  (32), 
and  neglecting  radiation,  yields 

%  =  T)]  (36) 

For  constant  collision  efficiency,  the  particie 
mass  equation  [Eq.  (33)]  will  integrate  immediately 
to  yield 

'■(O  ^  ro-^t  (37) 

where  rg  is  the  initial  particle  radius.  Equation  (36) 
is  a  linear,  first-order,  ordinary  differential  equation 
with  variable  coefficients  (since  r=  r(t})  which  can 
be  simply  integrated  using  an  integrating  factor  to 
obtain 


Fig.  1.  Comparison  of  analytical  and  numerical 
solutions  for  particle  temperature  history. 


The  analytical  and  numerical  solutions  are  virtu- 
aliy  indistinguishable  from  each  other,  except  at 
the  smaller  particle  sizes.  The  numerical  solutions 
were  performed  with  a  constant  (10“^  sec)  step 
size,  which  was  not  small  enough  to  capture  the 
fast  rise  for  the  smallest  particles.  This  result  vali¬ 
dates  the  bulk  of  the  numerical  solution.  The  only 
potential  errors  are  in  the  evaluation  of  the  non- 
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constant  properties  and  in  the  treatment  of  the 
radiation  term.  Note  in  Fig.  1  that  the  ultimate  tem¬ 
perature  of  the  particle  is  independent  of  the  initial 
particle  size.  The  initial  particle  size  affects  only  the 
rise  time. 

Including  the  radiative  term  in  Eq.  (36)  yields 


and  so  we  can  define 

T*=  lim[T{t)]  (43) 

t-^t* 

The  burnout  temperature  ( 1*)  exhibits  two  dis¬ 
tinct  behaviors  depending  upon  the  value  of  the 
exponent  5: 


where  A>is  now  treated  as  a  temperature-indepen- 
dent  constant.  Figure  2  presents  the  same  numeri¬ 
cal  results  as  Fig.  1,  with  the  addition  of  the  results 
from  the  numerical  solution,  including  radiation. 
The  analytical  solutions  have  been  removed  for 
clarity.  Note  that  neither  the  ultimate  particle  tem¬ 
perature  nor  the  particle  lifetime  change  with  the 
inclusion  of  the  radiation  term.  Inclusion  of  the  radi¬ 
ation  term  only  increases  the  time  required  to 
reach  the  ultimate  temperature. 

2590 

^  2580 

£  2570 

I  2560 

a  2550 

I  2540 

.S  2530 

1  2520 
<0 

“■  2510 

10-5  10-^  10-5  10-2  10-1 

Time,  sec 

Fig.  2.  Comparison  of  constant  property  numerical 
solutions  with  and  without  radiation. 

Critical  Behavior  Analysis 


T1  \ 


/ 

/ 


Initial  Particlejremp:  2500  K 

“  ‘  10~atrn  Oxygen  Pressure 
Roth  Oxidation  Rates 


lOP^  Initial  Particle  Diam,  nm 


Constant  Property  Numerical 
Solution  No  Radiation 
----  Constant  Property  Numerical 
Solution  With  Radiation 


The  presence  of  the  A^^^term  in  Eq.  (38)  indi¬ 
cates  a  singularity  at  burnout  (i.e.,  at  a  =  The 
burnout  temperature  f  is  the  particle  temperature 
as  this  singularity  is  approached.  The  time  at  which 
the  singularity  is  reached  {t=t)  can  be  found  by 
setting  r(t)  =  0\r\  Eq.  (37): 


Zfm 


(42) 


5< 

(44) 

5  ^  0'=^  T  *—>  oo 

(45) 

Allowing  the  collision  efficiency  (/)  to  vary  while 
considering  all  other  parameters  fixed,  we  can 
solve  Eq.  (40)  for  the  critical  value  of  f=  /that  cor¬ 
responds  to  5  =  This  yields 


f*  = 


(46) 


which  for  the  numerical  values  of  the  other  param¬ 
eters  yields  /  =  0.401.  This  means 


/</»=>  5  <  d7=>  (47) 


f>f*^^>0^T*-^oo  (48) 

Physically,  the  collision  efficiency  f  must  be 
between  0  and  1.  Therefore,  both  behaviors  are 
possible,  and  a  bifurcation  exists  at  f  =  f.  Results 
from  Eq.  (38)  are  plotted  in  Fig.  3  for  various  val¬ 
ues  of  the  collision  efficiency  f.  The  computations 
stop  one  time  step  before  burnout  to  avoid  the  infi¬ 
nite  temperature  at  the  singularity. 


Fig.  3.  Analytical  solution  for  temperature  his¬ 
tory  as  a  function  of  collision  efficiency. 
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Note  the  change  in  the  behavior  of  the  slope  of 
the  temperature-time  history  as  /becomes  larger 
than  f.  The  slope  at  burnout  changes  from  practi¬ 
cally  zero  to  very  large  as  the  collision  efficiency 
increases.  Figure  4  is  a  re-scaled  version  of  Fig.  3 
to  show  the  complete  curve  at  f=  0.41 .  The  “small 
f  behavior  would  be  typical  under  most  physical 
conditions.  The  “large  /”  behavior  would  be 
approached  only  under  the  conditions  of  very  effi¬ 
cient  reactive  collisions  (such  as  with  atomic  oxy¬ 
gen)  and  no  nonreactive  gases  (such  as  nitrogen). 
These  conditions  are  not  typically  found  in  flames, 
shock  tubes,  or  other  practical  combustion 
devices.  However,  these  conditions  may  be 
encountered  in  the  combustion  of  ablative  materi¬ 
als  (such  as  carbon)  in  atmospheric  reentry. 


Fig.  4.  Analytical  solution  for  temperature  history  as 
a  function  of  collision  efficiency  -  rescaled. 


The  large  temperatures  indicated  are  not  physi¬ 
cal,  since  other  processes  that  are  ignored  here 
would  become  important.  Two  such  processes  are 
thermal  vaporization  (which  will  continue  to  be 
ignored)  and  radiation.  Radiation  will  now  be  con¬ 
sidered. 

Equation  (41)  is  the  particle  energy  equation 
including  radiative  loss.  Figure  5  shows  the  numer¬ 
ical  solution  of  Eq.  (41)  for  particle  temperature 
versus  time.  The  collision  efficiency  (^  is  varied 
while  considering  all  other  parameters  fixed. 

The  “small  f”  and  “large  /”  behaviors  noted  in 
Fig.  3  are  also  seen  when  including  radiative  heat 
loss.  Figure  6  is  a  rescaled  version  of  Figure  5  to 
show  the  complete  curve  for  f=  0.41 .  Note  that  the 
addition  of  the  radiative  term  did  not  change  the 
ultimate  particle  temperature.  Once  again,  the  high 


Fig.  5.  Numerical  solution  for  temperature  history 
as  a  function  of  collision  efficiency  (with 
radiation). 


Fig.  6.  Numerical  solution  for  temperature  history 
as  a  function  of  collision  efficiency  (with 
radiation)  -  rescaled. 


temperatures  indicated  here  are  not  physical, 
since  thermal  vaporization  will  become  important 
at  about  4000  K.  This  mechanism  would  remove 
mass  (and  thus  energy)  from  the  particle  while 
holding  the  particle  temperature  down.  As  in  the 
previous  figures,  the  calculations  are  stopped  one 
time  step  before  burnout  to  avoid  the  singularity 
which  occurs  when  the  particle  radius  is  reduced 
to  zero. 

Figure  7  shows  the  solutions  for  both  the  ana¬ 
lytical  solution  without  radiation  and  the  numerical 
solution  with  radiation.  For  small  collision  efficien¬ 
cies,  the  curves  are  practically  indistinguishable, 
but  the  radiation  term  does  increase  the  tempera¬ 
ture  rise  time,  as  noted  in  the  previous  section.  For 
larger  values  of  the  collision  efficiency,  the  slope  of 
the  temperature  at  burnout  is  practically  zero  for 
the  no-radiation  cases,  but  is  very  large  for  the 
cases  considering  radiation.  The  inclusion  of  radia- 
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Fig.  7.  Comparison  of  analytical  and  numerical 
(with  radiation)  solution  for  temperature 
history  as  a  function  of  collision  efficiency. 


tion  reduces  the  value  of  the  collision  efficiency 
required  to  produce  a  large  slope  at  burnout.  The 
critical  value  for  infinite  temperature  at  burnout  is 
unaffected  by  the  inclusion  of  radiation. 

Summary/Conclusions 

A  bifurcation  is  seen  in  the  behavior  of  small 
particles  combusting  under  free  molecular  condi¬ 
tions.  For  relatively  low  collision  efficiencies  (or 
kinetic  rates),  the  particle  temperature  is  bounded 
as  the  particle  size  reduces  to  zero  at  burnout.  As 
the  collision  efficiency  is  increased,  the  particle 
temperature  at  burnout  increases,  but  still  remains 
bounded.  At  some  critical  value  of  collision  effi¬ 
ciency,  the  particle  temperature  at  burnout  is 
unbounded,  and  increases  without  limit.  Under  the 
assumption  of  constant  thermal  properties  of  both 
the  particle  and  the  surrounding  gas,  this  critical 
value  of  the  collision  efficiency  is  seen  to  be  a  sim¬ 
ple  function  of  the  gas  and  particle  specific  heats 
and  the  amount  of  mass  removed  per  collision. 
Calculations  were  presented  neglecting  and  includ¬ 
ing  the  effect  of  particle  radiation.  The  inclusion  of 
radiation  increases  the  time  required  to  reach  any 
given  particle  temperature,  but  does  not  substan¬ 
tially  change  the  critical  behavior;  i.e.,  the  particle 
temperature  at  burnout  becomes  unbounded  at  the 
same  critical  value  of  collision  efficiency,  whether 
radiation  is  included  or  neglected. 

The  impact  of  this  critical  behavior  will  be 
important  in  situations  involving  high  collision 
efficiency — i.e.  high  particle  temperatures  in 


extremely  oxidizing  conditions.  For  instance,  the 
use  of  laser  induced  incandescence  as  a  soot 
diagnostic  is  typically  applied  in  reducing  or  at 
most  weakly  oxidizing  conditions.  The  kinetic  time 
scales  are  long  enough  that  reactions  on  the  parti¬ 
cle  surface  during  the  laser  pulse  and  subsequent 
particle  cooldown  are  typically  Ignored.  If  soot  (or 
other  small  combustible  particles)  were  present  in 
an  oxidizing  atmosphere  (e.g.  ablation  or  inten¬ 
tional  particle  injection  into  the  boundary  layer  of  a 
re-entry  vehicle)  and  laser  induced  incandescence 
used  as  a  particle  concentration  diagnostic,  the 
potential  for  critical  vs.  non-critical  combustion 
behavior  would  have  to  be  assessed. 

For  the  carbon  particles  considered  here,  ther¬ 
mal  vaporization  (sublimation)  is  expected  to  have 
a  major  impact  on  the  critical  behavior  as  the  parti¬ 
cle  temperatures  approach  the  vaporization  tem¬ 
perature  of  about  4000  K.  Particle  vaporization  will 
be  added  to  the  model  in  the  near  future. 
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